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ABSTRACT 

We present: 1) a kinematic and morphological study of the giant Lya nebula associated 
with the radio galaxy MRC 2104—242 (z = 2.49) based on integral field spectroscopic 
VIMOS data from VLT ; 2) a photometric study of the host (proto?) galaxy based on 
Spitzer Space Telescope data. 

The galaxy appears to be embedded in a giant (> 120 kpc) gas reservoir that 
surrounds it completely. The kinematic properties of the nebula suggest that it is 
a rotating structure, which would imply a lower limit to the dynamical mass of ~ 
3 x 10 11 Mq. An alternate scenario is that the gas is infalling. Such a process would be 
able to initiate and sustain significant central starburst activity, although it is likely 
to contribute with less than 10% of the total stellar mass of the galaxy. 

The near- to mid-IR spectral energy distribution of the radio galaxy suggests the 
existence of a reddened, E(B — V) = 0.4 ±0.1, evolved stellar population of age > 1.8 
Gyr and mass (5 ±2) x 10 11 M Q . The implied formation redshift is Zf > 6. This stellar 
mass is similar to the stellar masses found for massive early-type galaxies at z ~2 in 
deep, near-infrared surveys. 
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1 INTRODUCTION 

In a study based on long-slit, medium resolution (R^700) 
Keck II/ VLT spectroscopy of a sample of 10 ultra steep 
spectrum, high redshift (z ~ 2.5) radio galaxies (HzRGs), 
we found that these objects are embedded in giant (often 
^100 kpc), low surface brightness halos of metal rich, ion- 
ized gas with quiescent kinematics (Villar-Martm et al. 2003 
[VM03]). By quiescent we mean that the kinematics do not 
appear to be perturbed by interactions between the radio 
and optical structures. The halos usually extend well be- 
yond the brighter, perturbed regions and sometimes beyond 



* Based on observations carried out at the European Southern 
Observatory, Chile (programs 073.B-0189, 065.P-0579 and 075. B- 
0729). and with the Spitzer Space Telescope. 



the radio structures. We proposed that the quiescent ha- 
los are the gas reservoirs from which the galaxies started to 
form and might still be forming. At some point, an active 
nucleus switches on and the halo becomes observable thanks 
to strong line emission powered by the ionizing continuum. 

Our work was seriously limited by the lack of spatial 
information in directions other than the radio axis, through 
which the long slits were aligned. For this reason, we are 
carrying out a programme of 3D integral field spectroscopy 
to characterize the morphological, kinematic and ionization 
properties in two spatial dimensions of the extended ionized 
gas in a sample of powerful radio galaxies at z ~2-3. 

In this letter, we present the results obtained for the 
radio galaxy MRC 2104-242 at z =2.49 (McCarthy et al. 
1990). The radio source has an angular extent of 24" (~200 
kpc). The WFPC2 and NICMOS HST images (Pentericci 
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et al. 2001) show that the host galaxy is very clumpy. A 
filamentary component more than 2" long is found aligned 
with the radio axis to within a few degrees. This region 
is likely to be undergoing a strong interaction between the 
radio and optical structures (VM03, Humphrey et al. 2005). 

The existence of an evolved stellar population in HzRGs 
has been subject of intensive study for the last two decades. 
The UV and optical restframe spectroscopic studies have 
only marginally detected the continuum, which may still 
contain significant non-stellar contribution (Vernet et al. 
2001). To detect the old stellar population we need to ob- 
serve at rest-frame near-IR wavelengths. This is one of the 
scientific goals of the Cycle 1 Spitzer survey of 70 high 
redshift radio galaxies (P.I. Stern, program ID 3329). We 
present results for MRC 2104—242 in this paper. 

A Q A = 0.73, Vim 0.27 and H 71 km s" 1 Mpc -1 
cosmology is adopted in this paper. 



2 OBSERVATIONS AND DATA REDUCTION 

The spectroscopic observations were made on UT 2004 June 
17 (ESO programme 073.B-0189(A)) using the Visible Mul- 
tiObject Spectrograph (VIMOS, Le Fevre et al. 2003), on 
the Nasmyth focus of the UT3 VLT unit. The instrument 
is equipped with an integral field unit made of 1600 mi- 
crolenses coupled to fibres, covering 27"x 27 /; on the sky with 
0.67" diameter fibres for the high spectral resolution mode. 
The HRbiue grating was used, with an effective wavelength 
range ^4150-6200 A and an instrumental profile 3.0±0.5 A. 
The total integration time was 4.5 hr (9 x 1800 sec). The see- 
ing FWHM during the observations was in the range ~1.0- 
1.4". 

The data were reduced using a modified version of P3d 
(Becker 2002), in addition to our own software optimized 
for VIMOS data (Sanchez & Cardiel 2005). The data were 
bias subtracted. The expected locations of the spectra were 
traced on a continuum-lamp exposure obtained before each 
target exposure. The fibre-to- fibre response at each wave- 
length was determined from a continuum-lamp exposure. 
Wavelength calibration was performed using arc lamp spec- 
tra and the telluric emission lines in the science data. After 
these basic reduction steps, a data cube was created for each 
exposure. The cubes were then recentred spatially by deter- 
mining the centroid of a nearby star at each wavelength. This 
recentering corrects for differential atmospheric refraction. 
The cubes were then combined using IRAF tasks, masking 
the broken and/or low sensitivity fibres. A 3a clipping algo- 
rithm removed cosmic rays. The sky background spectrum 
was estimated before subtraction by averaging spectra of 
object free areas of the data cube. The flux calibration was 
done using the long slit FORS-VLT spectrum discussed in 
VM03. Cross-talk effects were found to be negligible. These 
were evaluated using the field star in the VIMOS field of 
view. The fibre to fibre contamination was found to be <5% 
for adjacent spectra in the CCD, dropping to less than 0.1% 
for the 3rd adjacent spectra. 

To overlay the radio core position (Carilli et al. 1997) on 
the Lya image, we used the NICMOS HST image, where the 
radio core was assumed to coincide with the IR peak. This 
is probably not the case if the nucleus is heavily reddened. 
The NICMOS peak position in the Lya image was then 
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Figure 1. Lya image of MRC 2104—242 . Contours correspond to 
~5, 10, 15, 20, 25, 35, 45, 60, 70, 80, 90 and 100 % of the maximum 
contour value. Pixels marked in green and pink colours were used 
to extract the spectra in Fig. 2. Ionization cones with an opening 
angle of 90° and a vertex at the radio core position (red 'x') are 
shown (black lines). The spectra of two intercone regions (red 
squares) are also shown. Lya: is detected in both cases. 

calculated using a bright star that is present in both images. 
The la uncertainty in the Lya-radio registration is ~0.3". 

Details about the observations and data reduction for 
the Spitzer data will be presented by Seymour et al. (in 
prep.). Briefly, for MRC 2104-242, we obtained 3.6, 4.5, 
5.8 and 8.0 fim images on UT 2004 October 27 using the In- 
fraRed Array Camera (Fazio et al. 2004), and a 16 /xm image 
in "peak-up only mode" on UT 2004 October 21 with the In- 
fraRed Spectrograph (Houck et al. 2004). The IRAC "basic 
calibrated data" exposures were mosaiced us ing the M0PEX 
software package (http://ssc.spitzer.caltech.edu/postbcd/) 
with the pixels being resampled by a factor of 2. We used 
SEXTRACT0R (Bertin & Arnouts 1996) to measure the pho- 
tometry in a matched 1" aperture defined by the 3.6 /mi im- 
age. To obtain total magnitudes, we applied aperture cor- 
rections from in flight point spread functions (Lacy et al. 
2005). 



3 RESULTS 

A Lya image (Fig. 1) was created by adding all the frames 
in the spectral direction of the data cube across the Lya 
emission line spectral profile, covering the ^4225-4260 A 
spectral range. The large equivalent width of the line in all 
positions where it is detected ensures that the contamina- 
tion by continuum emission is negligible everywhere across 
the image. Fig. 1 shows two spatial components (SW and 
NE) separated by ^6 arcsec as previously described by Mc- 
Carthy et al. (1990). The radio core is located between the 
two clumps, where Lya presents a minimum. Similar Lya 
morphologies have been observed in other HzRGs (e.g. Reu- 
land et al. 2003). The SW component overlaps partially with 
the filamentary structure described in §1. 

Low surface brightness Lya emission is detected in re- 
gions that seem to be outside any plausible ionization cones. 
Two spectra were extracted from the two intercone regions 
shown in Fig.l, where we have assumed a typical cone open- 
ing angle of 90° and the position of the cone vertex coincid- 
ing with the radio core. Lya, although noisy, is clearly de- 
tected in the two intercone regions. Any rotation of the axis 
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Figure 2. 1-dimensional Lya: VIMOS spectra extracted from the 
NW and SE clumps, using the regions shown in Fig.l. The flux 
levels have been shifted for clarity. The clear split observed in the 
SW Lya profile is due to kinematics. We argue that the SWR 
and NE components are both quiescent gas. 



or reasonable shift of the vertex position leads to some of the 
high surface brightness emission lying outside the ionization 
cones. McCarthy et al. (1990) mentioned the detection of a 
diffuse Lya halo that surrounds the entire object. If the gas 
is ionized (but see Villar-Martm, Binette & Fosbury 1996), 
the ionization mechanism is an open question. Hot young 
stars and cooling radiation are two interesting possibilities. 



3.1 Lya kinematics 

Our goal is to map the 2D kinematic field of the quiescent 
gas (i.e. non perturbed by jet-induced shocks) and use it as 
a source of information about its origin and the formation 
process of the galaxy. If the gas motions are gravitational, we 
can constrain the dynamical mass of the system. Because of 
the high sensitivity of Lya to absorption by HI, it is often an 
unreliable kinematic tracer. Given the low S/N of the data, 
we cannot use other reliable tracers, such as the HeIIA1640 
line. Instead, we have compared the integral field results here 
with our earlier kinematic study (VM03) based on FORS2- 
VLT moderate resolution (6 A) long slit spectroscopy, which 
made use of both Lya and Hell. 

We show in Fig. 2 the 1-dimensional VIMOS spectra 
of the SW and NE clumps using two apertures (see Fig.l) 
matching as much as possible apertures 4 an d 3 of VM03. 
Lya in the SW clump is split into two components (blue 
(SWB) and red (SWR)), both of which are well repre- 
sented by Gaussian profiles. The line FWHM are 470±20 
and 550±20 km s _1 respectively. They are separated by 
820±20 km s _1 . The line in the NE component is well rep- 
resented by a Gaussian profile of FWHM 630±20 km s _1 
redshifted by 250±10 km s _1 relative to the SWR compo- 
nent. These results are in good agreement with those ob- 
tained with the FORS2-VLT data. In particular, the Hell 
line in the VM03 long-slit spectrum of the SW clump shows 
the same split as seen in the current VIMOS Lya profile. 
Since the Hell line is much less susceptible to absorption, 
this implies that the double structure is due to kinematics 
rather than absorption. 

We argue that the red Lya component from the SW 
clump (SWR) and the Lya emission from the NE clump are 
emitted by quiescent gas. The remarkable similarity in the 
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Figure 3. Overlay of the VIMOS Lya intensity contours of the 
quiescent gas (i.e., components NE and SWR) with the Lya ve- 
locity field (left) and line width a (right). The position of the 
radio core is indicated with a cross. Zero velocity corresponds to 
a redshift z = 2.490. The velocity field appears symmetric and 
ordered, suggesting either rotation or radial flows. 



kinematic properties (see below) of these two components 
strongly suggest that the gas in both regions follows a similar 
kinematic pattern with the same origin. On the other hand, 
the sharp and large velocity shift of the SWB component 
suggests that this is a different, decoupled kinematic com- 
ponent. Although detected almost all over the SW clump, 
this component is strongest relative to the SWR compo- 
nent in pixels adjacent and coincident with the filamentary 
structure (§1), where strong evidence for jet-gas interactions 
exist (Humphrey et al. 2005), suggesting that it is emitted 
by perturbed gas. 

We find that the Lya profile of the quiescent gas is 
not strongly distorted by absorption. VM03 found that in 
general, although the kinematic measurements (FWHM and 
velocity shifts) for Lya and Hell can be rather different for 
the perturbed gas (usually broad lines), they are are in good 
agreement within the errors for the quiescent gas. This is also 
the case for MRC 2104-242 (see fig.7, VM03). Therefore, 
the quiescent Lya emission is less affected by absorption, 
probably due to the narrowness of the lines. 

In summary, the Lya emission from the quiescent 
gas (components NE and SWR) can be isolated in 
MRC 2104—242 at each spatial location and it can be used 
as a reliable kinematic tracer. 

We now analyse the velocity field of the quiescent gas. 
The Lya profile was fit at each spatial pixel with Gaussian 
profiles. The results are shown in Fig. 3 (left) with the in- 
tensity contours of the quiescent gas overplotted. The errors 
are <30 km s _1 for most pixels. The velocity dispersion a 
map (corrected for instrumental broadening) is also shown 
(Fig. 3, right). The error is estimated to be <60 km s _1 . 

The kinematic pattern is rather symmetric and appar- 
ently ordered. The overall velocity field is very suggestive of 
a rotating structure, for which the half amplitude of the ro- 
tation curve is v ro t ^150 km s _1 . The a map is quite feature- 
less with a ^200-280 km s _1 for most spatial positions, with 
little noticeable difference between the two clumps other 
than a slight broadening of the line near the radio core. 
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If the active nucleus coincides with the dynamical center of 
the system, this line broadening further supports rotation. 



3.2 Spitzer data 

Table 2 presents the Spitzer photometry for MRC 2104-242, 
as well as near-infrared photometry from ground-based 
(J S ,K S ) ISAAC- VLT (ESO programme 075.B-0729(B)) and 
HST (F160W) images. We have corrected the near infrared 
photometry for emission line contamination using archival 
VLT/ISAAC JHK spectroscopy (ESO programme 65.P- 
0579(A)). We determine the contributions to be 15±5%, 
76±5%, and 51±5% to the total flux in J s , H and K s , re- 
spectively. 

Fig. 4 plots the near- to mid-IR spectral energy distri- 
bution (SED). The SED shows a clear rise between the K— 
and IRAC 3.6 /im bands, and again between 8.0 and 16 /mi. 
We now argue that the 3.6 to 5.8 /im points are dominated 
by an old stellar population. Figure 4 shows that the up- 
turn in the SED between 8.0 and 16 /im (see Table 2) is 
unlikely to be explained by a stellar population; the 16 /im 
point is most likely dominated by warm dust emission asso- 
ciated with the active nucleus. We assume the 16 /im point 
can be fit with a simple power law of slope —3 (Fig. 4), as 
determined from a sample of 12 radio galaxies with Spitzer 
16 and 24 /im data (Seymour et al., in prep.). Scaling this 
power-law to the 16 /im point implies that we can safely ig- 
nore any thermal dust emission at A G bs ^ 8.0 /im. This is 
further supported by the 8.0 /im upper limit. 

Similarly, the JHK points allow us to constrain the 
combined contributions due to scattered quasar light, neb- 
ular continuum and young stellar populations. These con- 
tributions all have a blue slope, while we detect a clear rise 
in the SED between 2.2 and 3.6 /im. We can therefore con- 
strain these combined contributions to <C5% in the IRAC 
bands. Note that the excess of the J band flux over the model 
predictions (Fig. 4) and over the H band measurement can 
only be explained by such contributions, since any old stellar 
population would be fainter in J than in H (especially when 
we invoke reddening) . Hence our strongest constraint comes 
from the H band detection. Based on restframe near-infrared 
spectroscopy of low-redshift active galaxies (e.g. Sosa-Brito 
et al. 2001), we expect emission line contributions to the 
IRAC bands to be <5%. 

We now use the H thru 8.0 /im SED of MRC 2104-242 
to derive the most probable stellar mass of the host galaxy 
(see, e.g., Yan et al. 2005 for a detailed discussion on un- 
certainties). The J-band point is omitted as it appears to 
be contaminated by light associated with either the buried 
AGN or a young stellar population, and the 16 /im point is 
ignored as its likely non-stellar in origin. We systematically 
compare the observed SED with a set of templates com- 
puted with PEGASE.2 models (Fioc & Rocca-Volmerange 
1997) via a x 2 minimization technique. Star formation his- 
tories (SFHs), ages, masses and dust extinction are free pa- 
rameters. We model simple stellar populations of different 
metallicities using a Salpeter (1955) IMF with a lower mass 
limit of 0.1 Mq and an upper mass limit of 12OM . We adopt 
SFHs that have been shown to consistently reproduce obser- 
vations of z ~ passively-evolving, early-type galaxies (Fioc 
& Rocca-Volmerange 1999). A more detailed description of 
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Table 1. Near to mid-IR photometry of MRC 2104-242 . The 
J S HK S photometry has been corrected for line emission (see 
text). 
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Figure 4. Near to mid-IR SED of MRC 2104-242 . The best 
fit is produced by a reddened, relatively evolved population of 
age 1.8 - 2.5 Gyr and stellar mass (5 ± 2) x 10 n M Q . The J s 
band point is dominated by a combination of scattered quasar, 
nebular continuum and young stellar emission. All of these con- 
tributions have a blue slope, and will therefore contribute <C5% 
to the IRAC photometry. A power law of index -3 (dotted line) 
represents the expected thermal dust emission. The figure shows 
that its contribution can be safely ignored below 8 /im. 



the fitting method will be presented in Rettura et al. (2005, 
in prep.). 

Results of the best-fit model are shown in Fig. 4. We 
estimate the errors for the mass by sampling the full proba- 
bility distribution in the multidimensional parameter space. 
We find a best-fit stellar mass estimate of (5 ± 2) x 10 n M© 
and a color excess E(B — V) = 0.4 ±0.1. This is similar 
to the stellar masses found for the most massive, early- type 
galaxies at z ~2 in deep, infrared surveys (e.g. Daddi et al. 
2004, Labbe et al. 2005). The best-fit age is 2.5 Gyr, which 
is very close to the age of the universe at the source redshift, 
(2.7 Gyr). The minimum age that still provides a reasonable 
(3cr) fit is 1.8 Gyr, implying a formation redshift Zf > 6. Fi- 
nally, we note that the SED cannot be fit by any simple or 
composite stellar population without invoking substantial 
reddening. 
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4 DISCUSSION AND CONCLUSIONS 

Let us assume that the giant Lya nebula in MRC 2104—242 
is supported by rotation. If its morphology maps the true 
gas distribution and the gas is settled on a disk, we derive 
an inclination angle i ~ 75 =b 5° relative to the plane of the 
sky and a dynamical mass ~ (3.2 ± 0.7) x 10 n M© within a 
radius of 7 arcsec (~ 60 kpc). This is similar to the stellar 
mass inferred from the Spitzer data (see §3.2), suggesting 
a surprisingly small dark matter fraction. However, the de- 
rived dynamical mass is likely to be a lower limit, since the 
gas is expected to be anisotropically illuminated by the hid- 
den quasar. Therefore, the disk inclination and the direction 
of the major axis cannot be determined (e.g., the gas we see 
might be close to the minor axis). 

If the giant halo in MRC 2104-242 is settled on to 
a well defined disk, then — — < 1, such that v ro t must be 

v rot 

> 240 km s -1 and M dyn >8x 10 n M o . On the other hand, 
it takes several rotation periods for gas to settle into a disk. 
If this is the case for the giant nebula in MRC 2104-242 
and we assume a minimum of three rotation periods (P ro t = 
27rr/i>rot) for a well-ordered disk, the 2.7 Gyr age of the 
universe at z = 2.49 requires v YO t > 400 km s _1 for which 
Mdyn > 2 x 1O 12 M . Smaller v rot would imply that the gas 
has not yet settled into a stable structure. 

It is possible that the SW and NE clumps are two dif- 
ferent objects rotating around a common center (e.g. De 
Breuck et al. 2005). However the fact that Lya emission sur- 
rounds the object completely suggests that it is a single, ro- 
tating structure. If the gas is rotating, this means that giant 
(r > 60 kpc) gaseous, rotating structures (disks?) can exist 
already at z ~ 2.5. These could be the progenitors of the gi- 
ant (> several tens of kpc) rotating disks discovered around 
several low redshift radio galaxies and early-type galaxies 
(e.g. Morganti et al. 2002). High redshift cooling flow models 
for galaxy formation (e.g. Haiman, Spaans & Quataert 2000) 
predict the formation of a giant, rotationally-supported disk 
as part of the process in less than 1 Gyr. Therefore, the 
evolved (> 1.8 Gyr old), stellar population revealed by the 
Spitzer data probably formed well before the disk. 

Alternatively, the general pattern of the velocity field 
is also suggestive of radial motions. The radio properties 
can discriminate between inflow and outflow. The NE radio 
lobe is brighter and more polarized (Pentericci et al. 2000). 
According to the Laing-Garrington effect (Laing 1988), this 
lobe is closer to the observer. Since the gas located at this 
side of the nucleus is the most redshifted, it must be infalling. 
The collapse of such a giant structure suggests that we are 
witnessing an early stage in the formation process of the 
central galaxy (e.g. Haiman, Spaans & Quataert 2000). 

Let us assume an infall velocity of ~ 150 x sin a km s _1 
and a radius ^ 60 x cos a kpc, where a is the angle between 
the ionized cone axis and the line of sight. Since this is a 
narrow-line radio galaxy, a is expected to be large, a > 
60°. Therefore, after < 2.5 x 10 8 yr, the 60 kpc halo should 
collapse and disappear in a very small redshift interval (by 
z ~ 2.2). The existence of giant, quiescent halos in radio 
galaxies at very different redshift s (van Ojik et al. 1996, 
Villar-Martm et al. 2005) suggests that they are long-lived 
structures, unless there are multiple episodes of Lya halos. 
The implication is that much larger reservoirs of invisible 



gas must exist and/or some mechanism must deccelerate 
the infall process at some point (Iono et al. 2004). 

As in VM03 and using the geometrical information pro- 
vided by the VIMOS Lya image, we estimate an ionized gas 
density n ~ 55 cm -3 (see VM03 for uncertainties on this cal- 
culation) and a gaseous mass (neutral and ionized) within 
60 kpc of - 1O 1O M . A mass deposition rate in the range 
~ (40 — 250) M© yr -1 is inferred. This infall rate might be 
able to initiate and sustain central starburst activity at sev- 
eral tens of Mq yr -1 (Iono et al. 2004). Since the object con- 
tains a ^ 1.8 Gyr old population as well as the young popu- 
lation triggered by the current inflow, at least two episodes 
of star formation are inferred for the formation of this mas- 
sive galaxy. If the halo has ^10 10 Mq, this is probably the 
maximum it will contribute to build up the stellar mass of 
the galaxy, so that the infalling gas will only be <10% of 
the galaxy mass. 

In summary, MRC 2104—242 contains an evolved stel- 
lar population that formed at z > 6. The kinematic proper- 
ties of the extended gas suggest that it is associated with a 
giant rotating structure, the remnant of the formation pro- 
cess of the galaxy. Alternatively, the gas might be infalling 
during an early phase in the galaxy formation process. Ex- 
tending a similar study to a larger sample of objects will 
be very useful to understand the demographics and nature 
of ordered kinematics in the quiescent gas associated with 
HzRGs. 
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